Simple Ginzburg-Landau Theory for Vortices in a Crystal Lattice 
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We study the Ginzburg-Landau model with a nonlocal quartic term as a simple phenomenological 
model for superconductors in the presence of coupling between the vortex lattice and the underlying 
crystal lattice. In mean-field theory, our model is consistent with a general oblique vortex lattice 
ranging from a triangular lattice to a square lattice. This simple formulation enables us to study 
the effect of thermal fluctuations in the vortex liquid regime. We calculate the structure factor of 
the vortex liquid nonperturbatively and find Bragg-like peaks with four-fold symmetry appearing in 
the structure factor even though there is only a short-range crystalline order. 

PACS numbers: 74.20.De, 74.60. Ge 
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It is of great interest to study vortex lattice structure 
and correlations in superconductors in the presence of 
coupling between the vortex lattice and the underlying 
crystal lattice. Various experimental probes including 
neutron diffraction |Q, Bitter decoration 0], and scan- 
ning tunneling microscopy have been used to reveal 
a range of vortex lattice structures from the usual trian- 
gular lattice to a general oblique lattice and a square lat- 
tice oriented along a specific direction of the crystal axis. 
An important feature of these structures is the emergence 
of the four-fold symmetry representing the symmetry of 
the underlying crystal lattice. This effect of the crys- 
tal lattice on the vortex lattice is found to be dependent 
upon the external field in such a way that a triangular 
vortex lattice is observed at low fields, while at higher 
fields a square lattice is observed Q] . 

In order to study the vortex lattice structure, one usu- 
ally uses a Ginzburg-Landau (GL) phenomenological the- 
ory. Since the usual GL theory is rotationally invariant, 
one needs additional terms that break this symmetry to 
account for the appearance of the four-fold symmetric 
vortex lattice structure. The conventional way to in- 
clude this effect is to introduce terms quadratic in the 
order parameter in the GL free energy with fourth or- 
der derivatives The observed four- fold symmetric 
vortex lattice structure can be explained within these 
formalisms. However, the equations involved in these 
theories are very difficult to handle even at the linearized 
level where one has to resort to approximate or numerical 
methods. 

In this paper we propose a much simpler phenomeno- 
logical model for vortices in a crystal lattice. Our model 
is the usual GL theory for a one-component complex or- 
der parameter, except that the term quartic in the order 
parameter, \P(r) is nonlocal. We consider the situation 
where the effect of the underlying crystal on the order 
parameter symmetry can be summarized into an appro- 
priate form for the nonlocal interaction potential. The 
vortex lattice structure is determined within our model 
without much calculational effort. We find that the four- 



fold symmetric vortex lattice structure can be modelled 
with the appropriate choice of a minimal number of pa- 
rameters describing the nonlocal quartic interaction. The 
a-b plane anisotropy found in high-temperature super- 
conductors is also incorporated in our model in the usual 
way through second order gradient terms. We believe 
that this model captures the same physics as in the higher 
derivative approaches of Ref. 0-0), but is just much sim- 
pler to handle calculationally. 

One feature of our simple formulation is that it allows 
a study of the effect of thermal fluctuations in the vor- 
tex liquid regime. This would be quite impossible in the 
conventional formulations. We apply the nonperturba- 
tive method developed by us in Ref. to calculate the 
structure factor of the vortex liquid, which is measured 
in neutron scattering experiments. One of the main re- 
sults of the present work is that one can observe in the 
structure factor the ring patterns expected for a liquid 
state but broken up into Bragg-like peaks. This sug- 
gests that even in the vortex liquid state where there is a 
short-range crystalline order, a weak four- fold symmetric 
coupling to the underlying crystal may produce the spots 
observed in neutron scattering experiments and usually 
attributed to the formation of the vortex crystal state. 

The model we study in this paper is based on the 
GL free energy with a nonlocal quartic interaction for 
a two-dimensional superconductor in a magnetic field 
B = V x A: 
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where a, (3 are phenomenological parameters, m x the ef- 
fective mass in the x-direction, m y in the y-direction, 
and D = — iV — (e* /hc)A. The a-b plane anisotropy 
is represented by the ratio of the effective masses. For 
later use, we define a = (mj/mj) 1 ' 4 . The only differ- 
ence between our model and the usual GL theory is the 
nonlocal quartic interaction term represented here by a 
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general function g(r) = g(x,y), which is equal to the 
delta-function S 1 - 2 ^ (r) for the usual local GL theory. 

One reason we study the two-dimensional version of 
the model in this paper is that one can easily apply the 
nonperturbative method in Ref. j?]] to calculate the vor- 
tex liquid structure factor. But, more importantly, as 
noted by one of us |s) , the phase correlation length par- 
allel to the field direction in a bulk superconductor grows 
exponentially as one approaches zero temperature. When 
this length scale becomes comparable to or larger than 
the sample size in the low temperature regime, the sys- 
tem effectively behaves as a two-dimensional thin film 
with phase coherence across the sample and a very low 
effective temperature We note that most of the ex- 
periments mentioned earlier are performed in this regime, 
and we expect that they can be described by the present 
model when the parameter a is set to large negative val- 
ues, i.e. low temperatures. 

A nonlocal quartic interaction as in (|l|) appeared in 
the renormalization group study || of this system near 
its upper critical dimension. Even if one starts from a 
local theory, renormalization always drives the quartic 
term into an effective nonlocal one. It is not our aim 
here to derive an explicit form of g(r) from a microscopic 
theory. Instead, we take (|l|) as our starting point for 
the phenomenological description of superconductors in 
the presence of an interaction between the vortex lat- 
tice and the crystal lattice, and show that the variety of 
vortex lattice structures observed in experiments can be 
explained using a very simple form of g(r). 

The main ingredient one has to incorporate into the 
phenomenological construction of g(r) is the presumed 
four-fold symmetry of the underlying crystal. (General- 
izations to other crystal symmetries are of course possi- 
ble). In the present work, we construct g(r) such that it 
contains a term with explicit four-fold symmetry in ad- 
dition to a rotationally symmetric term. Thus, we take 
for the Fourier transform g(k) = J d 2 rg(r) exp(ik • r), 

g(k) = exp [ - C(-) 4 (l - e cos4(0 + 0„))] , (2) 



where k = (k cos 9 , k sin 8) and /i = */ e*B/hc is 
the inverse magnetic length. In particular, g(k) = 
exp[-C{(l - e){ki + fc 4 ) + 2(1 + 3s)k^}/ f x i ] for 8 = 0. 
In (0), we introduced three parameters C, e and 8q in 
such a way that the overall constant C controls the 
strength of the interaction between the vortex lattice and 
the underlying crystal (note that when C — ► 0, a local 
GL theory is recovered), and the dimensionless parame- 
ter e (0 < e < 1) represents the strength of the four- fold 
symmetric interaction compared to the rotationally sym- 
metric one. As will be discussed later, O controls the 
orientation of the vortex lattice with respect to the un- 
derlying crystal. This choice for g is certainly not unique. 
One could use various other forms with a four-fold sym- 



metric term as long as they are positive definite (to en- 
sure the stability of GL theory) and they do not introduce 
any unphysical non-analyticity (Note that in order to en- 
sure the analyticity of g(k) near k = with the cos AB 
term, we need at least fc 4 terms in (^)). But, in the pres- 
ence of the growing length scale mentioned earlier, we 
might expect that universality applies and the present 
form of the nonlocal quartic interaction does not alter 
the physical results. Note that we measure the wave- 
vectors with respect to the inverse magnetic length. But 
this is just for convenience in later calculations. A more 
natural length scale for the nonlocal kernel will be the 
lattice spacing of the crystal lattice, Iq. Therefore, the 
dimensionless parameter C which appears in (|^) will de- 
pend on the ratio of two length scales, i.e. C ~ (Iqh) 4 - 
As the magnetic field increases, the inter- vortex spacing, 
fi , gets smaller, so the nonlocal interaction term be- 
comes more important. This qualitative feature of our 
model is consistent with experimental findings where the 
four-fold symmetric vortex lattice structure is observed 
only in the high field regime. 

It is convenient to map the free energy functional of 
Eq. (Q) into a form with isotropic gradient terms using 
the following transformations: a(r') = (a x (r'),a y (r')) = 
(<j~ 1 A x (r), aA y (r)), and h(r') = <?(r), where r' = 
(x',y') = (ax,a~ 1 y), or the Fourier transform h(k') = 
lj{<jk x ,(j~ x k!^). Then Eq. (||) becomes in the new order 
parameter ip(r') = ^(r) = ty(cr~ l x' , ay') : 

F[i>] = J dV(|l|DV| 2 + a|^(r')| 2 ) 
+ dMdV^WrMri-r^l^)! 2 , (3) 

where D' = — iV' — (e* /hc)a(r') and m = [m x m y ) 1 l' 1 . 

Within mean field theory, the structure of the vortex 
lattice in our model can be determined in exactly the 
same way as in Abrikosov's work using the lowest Lan- 
dau level approximation (LLL) M|. We look for a pe- 
riodic solution to the linearized GL equations while re- 
stricting the order parameter to the space spanned by the 
LLL wavefunctions. In the Landau gauge, a = (—By' , 0), 
the normalized solution quasi-periodic over the two pe- 
riodicity vectors, r[ = 1(1,0), r' u = l(C,r]) is given by 
V(r') ~ cf>(r'\0) ee (2r 1 )^exp(-^y'y2)8 3 (TT(x' + 
iy')/l\C + iv) with the theta function 83. The magnetic 
length is fixed by the flux quantization condition; 
27t/j~ 2 =(area of unit cell)= l 2 r\. A useful representation 
for <p(r') is @ 

OO 

|0(r')| 2 = J2 (-ir"exp(-G 2 /V + zG.r'), (4) 

m.n— — 00 

where G = /j, 2 l(r]m, n — Cm) is the reciprocal lattice vec- 
tor corresponding to the periodicity vectors in the r'- 
space. The mean- field free energy density is given by 



2 



■"MF 



-ol 2 h /2(3(3a, where o>h = a + Tie*B/2mc 



determines the upper critical field H C 2(T), and the 
Abrikosov ratio (3a, which accounts for the contribution 
from the quartic term, is given by 



Pa = ]Tcxp(-G72 M 2 ) h(G). 



(5) 



In order to find a minimum free energy configuration, one 
has to minimize Pa- 

Since the sum in (|5|) converges very quickly, it is not 
difficult to find a configuration that minimizes Pa for 
given values of C, e and 9$. In the limit where C — ► (the 
low-magnetic field regime), one recovers the local GL the- 
ory, and Pa attains the familiar minimum value, 1.159 • • • 
for a triangular lattice in the r'-space. Also, since rota- 
tional invariance exists in the r'-space when C — ► 0, the 
triangular lattice has no preferred orientation with re- 
spect to the underlying crystal. When transformed back 
to the original r-space, it results in a distorted triangular 
lattice due to the a-b plane anisotropy. 

In the high field regime where one cannot neglect the 
effect of the nonlocal quartic term, the situation is differ- 
ent. Because of the a-b plane anisotropy carried over to 
/i(k'), the system in the r'-space is neither rotationally 
invariant nor four-fold symmetric. Therefore, the vortex 
lattice is orientated with respect to the crystal. In our 
model, the parameter 8q specifies the orientation. In gen- 
eral, we obtain an oblique lattice whose form depends on 
the values of C and e. In the original r-space, this lattice 
is further distorted due to the a-b plane anisotropy. To 
be more specific, we parametrize the periodicity vectors, 
r[ and rL, in terms of a centered rectangular lattice for 
which one can write £ = 1/2 and rj — (1/2) tan^', where 
<p' is an angle between r{ and r' u . We consider a general 
case where this lattice is rotated by an angle tp' . The re- 
sulting lattice can also be regarded as an oblique lattice 
with two primitive vectors of equal length and an angle 
2ip' (or 7r — 2(p') between them. For simplicity we focus 
on the case where 9q = 0. For given C, e and 9q = 0, 
we look for <p' and ip that minimizes Pa- For 8q = 0, 
we find that the minimum free energy configuration al- 
ways corresponds to <p' = 0, where one of the primitive 
vectors, rj coincides with the x'-axis. Any other orienta- 
tion of the vortex lattice can be obtained using different 
values of 9q- The angle tp' that gives the minimum free 
energy changes continuously from ~ 60° corresponding 
to a distorted triangular lattice to ~ 45° for a distorted 
square lattice as the four-fold symmetric coupling e in- 
creases from toward 1 for fixed C. We find that for C 
greater than some value C c ~ 0.015, there exists e c which 
depends on C such that the vortex lattice remains as a 
distorted square lattice for e > e c . A similar behavior 
to this was obtained in Ref. [Ej. Now, since the original 
order parameter ^(r) is quasi-periodic with respect to 
Z(er _1 ,0) and Z(<7 -1 C, C7?) for <p' — 0, the oblique lattice 



in the r-space has an angle ip, where tan <p — o~ 2 tamp' . 
(This relation becomes more complicated if p' ^ 0.) To 
summarise, within mean-field theory, the structure of the 
vortex lattice is mainly determined by the anisotropy pa- 
rameter a and the four-fold symmetric coupling e, and 
the orientation by 0$. All these parameters can in prin- 
ciple be fixed by experiments which determine the flux 
lattice structure as a function of the magnetic field. 

Thermal fluctuations, which are especially important 
in high-T c materials, melt the mean-field vortex lattice 
into a vortex liquid. For the local theory (C = 0), 
the effect of thermal fluctuations around the mean-field 
solution was studied by Eilenberger jll| using the or- 
thonormal basis for the LLL wave functions, <P{y\vq) = 
exp(i/i 2 iE2/o) < H r |0), where r spans one fundamental cell, 
or q = /i 2 (yo, ~xo) belongs to the first Brillouin zone. 
There are two different modes of excitation, whose en- 
ergies are denoted by e±(q). In the long wave length 
limit, q — ► 0, the hard mode behaves as e+(q) — 
const + 0(q 2 ) while the soft mode takes the form e-(q) = 
(a /2)(<7 4 //i 4 ) + 0(q 6 ) as q -> 0. The soft mode cor- 
responds to an incompressible shear deformation of the 
vortex lattice Jl3| and ao can be identified with the shear 
modulus C66 of the triangular lattice. 

When terms which break the rotational symmetry are 
present, the orientation of the vortex lattice is locked to 
the underlying crystal. Therefore we expect that there 
exists an excitation energy cost associated with a rigid 
rotation of the vortex lattice against the crystal lattice. 
We shall calculate this energy when the coupling to the 
underlying lattice is very weak, i.e. C < 1. Follow- 
ing Ref. pH , we first determine the soft mode energy 
e_(q) associated with ([!]) to the lowest order in C. (We 
assume a = 1 for simplicity.) After lengthy but oth- 
erwise straightforward algebra, we obtain the following 
anisotropic expression: 



e-(q) 
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(1 + C ai )(q 4 /^) 



+ Ce(P iq i + 2p 2 q 2 x q 2 y + P 3 q 4 y )/^\ + 0(q 6 ,C 2 ) (6) 

with known numerical constants a.\, Pi, ffa, and P3. This 
corresponds to a general form of the effective free energy 
p4[| for the displacement u, which should involve the ro- 
tation fields, Vij = [piUj — djUi)/2 as well as the usual 
strain fields, Uy = [diUj + djUi)/2. In our case where the 
vortex lattice is incompressible and four-fold symmetric, 
the effective free energy density reduces to 

-{\i{d x u x - dyu y ) 2 + 4:X 2 u 2 xy + 4uw^ + 8£u xy v xy } 

with four energy constants, Ai,A2,w, and £. For C = 0, 
one only has the purely elastic part (u> = £ = 0) and 
Ai = A2 = C66- For small C, the shear modulus will 
have a O(C) correction, and the rotation modulus u) and 
the coupling £ between the rotation and strain fields will 
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j), we obtain using the 

- Cat + Ce{2(3 2 - ft - 
c 66 Ce(/3i + ft)/2, and 



be proportional to Ce. From ( 
method of Ref. [jl3| Ai = cgq(1 
ft)/4), A 2 = c 66 (l + Con), uj = 

£ = cfcaCeGSi - #0/2. 

A useful quantity in studying the effect of thermal fluc- 
tuations is the so-called structure factor A(k) of the vor- 
tex liquid, which is proportional to the Fourier transform 
of the density-density correlation function (|5 , (r)| 2 |5'(r + 
R)| 2 ) c . In the low temperature regime, the length scale 
l± governing the degree of a short-range crystalline order 
perpendicular to the magnetic field in the vortex liquid 
becomes very large (l± 3> Iq). We expect then that a tiny 
amount of coupling to the underlying lattice might be 
able to break the rotational symmetry of the vortex liq- 
uid system even though there is only a short-range trans- 
lational order. For example, a ring-like pattern expected 
in a structure factor will be broken up into Bragg-like 
spots even in a liquid state (these spots will not be delta- 
function peaks, i.e. not true Bragg peaks). This is to be 
contrasted with the usual explanation of the appearance 
of Bragg-like spots in neutron diffraction patterns using 
a phase transition from a vortex liquid state to a vortex 
lattice state. We can estimate the angular dispersion 86 
of these spots using the above discussion on the rotation 
modulus: Since the energy associated with a rigid rota- 
tion by 89 of the crystalline region of area l\ is given by 
l"^_ui(S8) 2 , by equating this to k B T, one finds that 



C and e that correspond to a moderately strong cou- 
pling between the vortex lattice and the underlying crys- 
tal (e = 0.5). One can clearly observe four bright spots, 
which correspond to the nearest peaks in the structure 
factor, emerging from a ring-like pattern. To observe the 
next and higher order spots (which will give a structure 
factor a closer resemblance to that expected for a tri- 
angular lattice) one has to obtain the structure factor 
at lower temperatures. We find in general that, as the 
coupling between the vortex lattice and the underlying 
crystal gets weaker, one has to go to lower temperatures 
to observe the Bragg-like spots. This fact is in qualita- 
tive agreement with (^) (although the present numerical 
calculation is not done in the strict weak coupling limit). 



(<56>) 2 ~ k B T/tull = {0M2)k B T/c 66 l 2 j_Ce, 



(7) 




where we have used the numerical values for ft and ft. 

As the temperature is raised, according to (Q), one 
needs larger values of Ce to observe the Bragg-like 
spots. At moderately low temperatures, an approxima- 
tion scheme, called the parquet resummation method |t]] 
is accessible for the calculation of the structure factor of 
the vortex liquid. Using this method, one can explic- 
itly observe a ring-like pattern in the structure factor 
is broken up into Bragg-like spots as the temperature 
is lowered. For an isotropic system described by (^), 
it is straightfoward to apply the parquet resummation 
method to calculate the structure factor A'(k') in the k'- 
space. The only difference compared to the usual local 
theory is that one starts from the bare quartic potential 
h(k') which explicitly breaks the rotational symmetry of 
the local theory. For the detailed form of the nonper- 
turbative equations one has to solve numerically for A', 
and the reader is referred to Ref. |t]] for details. The 
structure factor in the original space is given simply by 
A(k) = A'(k'), where (k' x ,k' y ) = {a- l k x ,ak v ). 

We have calculated A(k) for various values of the pa- 
rameters, C, e and at different temperatures down to 
olt ~ —6.7, where the temperature is represented by the 
dimensionless quantity, ar = ct£r\/27r//3^ 2 , which goes 
to — oo as one approaches zero temperature. Fig. ^ shows 
a contour plot of A(k) at ar — —6.3 for the values of 



FIG. 1. A contour plot of A(k) at ar — —6.3. The param- 



eters used are C = 0.01, e = 0.5, do = 0, and a 



1.15. 



In summary, we considered a simple phenomenological 
model for vortices in a crystal lattice using a nonlocal GL 
theory. As well as explaining the observed four-fold sym- 
metric vortex lattice structures within mean-field the- 
ory, the present model suggests that there is a possibil- 
ity of observing Bragg-like spots within the vortex liquid 
regime as a consequence of coupling of the vortices to the 
underlying crystal. 
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